ASYMPTOTICS OF MULTIVARIATE SEQUENCES II. 
MULTIPLE POINTS OF THE SINGULAR VARIETY. 
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f5 I Abstract. Let -F'(z) — '^^ Orz"" be a multivariate generating function which is meromorphic in 

^ ■ some neighborhood of the origin of C*, and let V be its set of singularities. Effective asymptotic 

expansions for the coefficients can be obtained by complex contour integration near points of V. 

In the first article in this series, we treated the case of smooth points of V. In this article 

we deal with multiple points of V. Our results show that the central limit (Ornstein-Zernike) 

behavior typical of the smooth case does not hold in the multiple point case. For example, when 

V has a multiple point singularity at (1, . . . , 1), rather than ttr decaying as |r|~^''^ as [r| -^ oo, 

\^ ' Or is very nearly polynomial in a cone of directions. 

■5 



O 



1. Introduction 

In |P W02bj ■ we began a series of articles addressing the general problem of computing as- 

K^ I ymptotic expansions for a multivariate sequence whose generating function is known. Such 

^^ ■ problems are encountered frequently in combinatorics and probability; see for instance Ex- 

f— ^ . amples 2 - 8 in Section 1 of |Pemj . which collects examples from various sources includ- 

VO I ing |LL99l IFM771 IWil94l ICom74j . Our aim is to present methods which are as general as 

O ■ possible, and lead to effective computation. Our apparatus may be applied to any function 

whose dominant singularities are poles. Among other things, we showed in |PWn2bj that for 

all nonnegative bivariate sequences of this type, our method is applicable (further work may be 

required in one degenerate case). 

The article |PWn2bj handled the case when the pole variety was smooth. The present work 
is motivated by a collection of applications where the pole variety is composed of intersecting 
branches. Several examples are as follows. 



o 









(1 - 2 x){l- 2y) 
arises in Markov modeling |Karf)2j : 

(2) 



xv 

2+ ^ 



1 - xy{l + x + y + 2xy) 
G(z) 



uUh^) 



F(z) 

where the Lj are affine, arises in queueing theory |BM93j : 
(3) 



x"-y''z'' 
F{x,y,z) 



(1 - x)il - y){l - z){l - xz){l - yz) 
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is a typical generating function arising in enumeration of integer solutions to unimodular 
linear equations, and is the running example used in |dLS03j : 

(4) 

^^^'^^ = (l-(l/3)x-(2/%)(l-(2/3)x-(l/3)y) 
counts winning plays in a dice game (see Example II .31) . 

(5) 

■z/2 

F(x,y,z = -— 

{l-yz)P{x,y,z) 

arises in the analysis of random lattice tilings |('EP96j . In this case, in addition to the 
factorization of the denominator, there is an isolated singularity in the variety defined by 
the polynomial P, for which reason the analysis is not carried out in the present paper 
but in tCP03; . 

These examples serve as our motivation to undertake a categorical examination of generating 
functions with self-intersections of the pole variety. As in the above cases, one may imagine 
that this comes about from a factorization of the generating function, although we show in 
Example 13.21 that the method works as well for irreducible functions with the same "multiple 
point" local geometry. 

Our approach is analytic. For simplicity, we restrict to the two- variable case in this intro- 
duction, though our methods work for any number of variables. Given a sequence ars indexed 
by the 2-dimensional nonnegative integer lattice, we seek asymptotics as r, s — > oo. Form the 
generating function F{z, w) = J2r s o-rsz'^w^ of the sequence; we assume that F is analytic in 
some neighborhood of the origin. 



The iterated Cauchy integral formula yields 

1 f f F{z,w) 



dw dz, 



where C and C are circles centered at and F is analytic on a polydisk containing the torus 
C X C' . Expand the torus, by expanding (say) C, slightly beyond a minimal singularity of F 
(that is, a point {zo,wo) at which the expanding torus first touches the singular set V of F). 
The difference between the corresponding inner integrals is then computed via residue theory. 
Thus Urs is represented as a sum of a residue and an integral on a large torus. One hopes that 
the residue term is dominant, and gives a good approximation to ars- The residue term is itself 
an integral (since the residue is taken in the inner integral only). A stationary phase analysis 
of the residue integral shows when our hope is realized, and yields an asymptotic expansion for 
a^s' 

In |PW02bj we considered the case when the minimal singularity in question is a smooth 
point of V. The present article deals with the case where the minimal singularity is a multiple 
point: locally, the singular set is a union of finitely many graphs of analytic functions. This 
case includes the smooth point analysis of |P W02bj . and Theorem 13 . 91 with n set to essentially 
generalizes the analysis of that article). In the two- variable case, we know from Lemma 6.1 of 
|PWn2bj that every minimal singularity must have this form or else be a cusp (see also |Tsi93[ 
Lemma 3.1]), though more complicated singularities may arise in higher dimensions. 

The rest of this article is organized as follows. In the remainder of this section we describe 
in more detail the program begun in |PWn2bj and continued here and in future articles in this 
series. Section |21 deals with notation and preliminaries required for the statement of our main 
results. Those results, along with illustrative examples, are listed in Section IHI Proofs are given 
in Section We discuss some further details and outline future work in Sectional 
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Details of the Program. Our notation is similar to that in |PW02bj . For clarity we shall 
reserve the names of several objects throughout. We use boldface to indicate a (row or column) 
vector. The number of variables will be denoted d+1. The usual multi-index notation is in use: 
z denotes a vector (zi, . . . , z^+i)'^ G C^"*"^, and if r is an integer vector then z"" = Y[j ^/ ■ We 
also use the convention that a function, ostensibly of 1 variable, applied to an element of C"'^^ 
acts on each coordinate separately — for example e^ = (e^^, . . . ,e^'''+i). Throughout, G and H 
denote functions analytic in some polydisk about and F = G/H = ^j.Orz'". The set where 
H vanishes will be called the singular variety of F and denoted by V. Let D(z),T(z) denote 
respectively the polydisk and torus (both centered at the origin) on which z lies. 

A crude preliminary step in approximating a^ is to determine its exponential rate; in other 
words, to estimate log|ar| up to a factor of 1 + o(l). Let V denote the (open) domain of 
convergence of F and let logP denote the logarithmic domain in M "'"^, that is, the set of x G M "''^ 
such that e^ G P. If z* G D then Cauchy's integral formula 



2m J Jt{z*) z''^^ 
shows that Or = Odz*!"""). Letting z* -^ dT> gives 

log |ar| < — r • log |z*| + o(|r|), 
and optimizing in z* gives log |ar| < 7(r) + o(|r|) where 

(1.2) 7(r) := — sup r-x. 

xSlogD 

The cases in which the most is known about a^ are those in which this upper bound is correct, 
that is, log |ar| = 7(r) + o(|r|). To explain this, note first that the supremum in (|1.2|) is equal to 
r • x for some x G d\ogD. The torus T(e'') must contain some minimal singularity z* G V n dT>. 
Asking that log |ar| ~ — r • z* is then precisely the same as requiring the Cauchy integral (|1.1|1 
— or the residue integral mentioned above — to be of roughly the same order as its integrand. 
This is the situation in which it easiest to estimate the integral. 

Our program may now be summarized as follows. Associated to each minimal singularity z* 
is a cone K(z*) C (]R+)'^+^. Given r, we find one or more z* = z*(r) G V Pi dT) where the 
upper bound is least. We then attempt to compute a residue integral there. This works only if 
r G K(z*) and if the residue computation is of a type we can handle. Our program is guaranteed 
to succeed in some cases, and conjectured to succeed in others. It is known to fail only in some 
cases where the a^ are not nonnegative reals (not the most important case in combinatorial or 
probabilistic applications) and even then a variant seems to work. 

To amplify on this, define a point z* G V to be minimal if z* G dD and each coordinate of z* 
is nonzero. There are only three possible types of minimal singularities |PW02bl Lemma 6.1]), 
namely smooth points of V, multiple points and cone points (all defined below). It is conjectured 
that for all three types of points, and any r G K(z*), we indeed have 

log lorl = 7(r) + o(|r|) = -r • log |z*| + o(|r|) . 

This is proved for smooth points in |PWn2bj via residue integration, and the complete asymptotic 
series obtained. It is proved in the present work for multiple points under various assumptions; 
the fact that these do not cover all cases seems due more to taxonomical problems rather than 
the inapplicability of the method. The problem remains open for cone points, along with the 
problem of computing asymptotics. 
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Figure 1 . Local picture of a double pole 

It is also shown in |PW02bj that when Or are all nonnegative, then r G K(z*(r)), and therefore 
that a resolution of the above problem yields a complete analysis of nonnegative sequences with 
pole singularities. 

When the hypothesis that Cr be real and nonnegative is removed, it is not always true that 
r G K(z*(r)). In all examples we have worked, there have been points z* ^ D for which 
r G K(z*) and a residue integral near z* may be proved a good approximation to Ur', the method 
is to contract the torus of integration to the origin in some way other than simply iterating 
the contraction of each coordinate circle to a point. Thus a second open question is to settle 
whether there is always such a point z* in the case of mixed signs (see the discussion in ;PW02Bl 
Section 7]). 

The scope of the present article is as follows. We define multiple points and carry out the 
residue integral arising near a strictly minimal, multiple point. Unlike the case for smooth points, 
the integral is not readily recognizable as a standard multivariate Fourier-Laplace integral, and 
a key result of this article is a more manageable representation of the residue to be integrated 
(Corollarv l4.Hl and Lemma l4.4|) . The Fourier-Laplace integrals arising fall just outside the scope 
of the standard references, and we have been led to develop generalizations of known results. 
These results, which would take too much space here, are included in |P Wn2aj . 

The asymptotics arising from multiple point singularities differ substantially from asymptotics 
in the smooth case. In the remainder of this introduction, we give examples to illustrate this. 

Example 1.1 (simplest possible multiple point). Let F{z,w) be a two-variable generating func- 
tion and suppose that the point (1, 1) is a double pole of F (thus F = G/H with G{z,w) ^ 
and H[z, w) vanishing to order 2 at (1, 1)). If F has no other poles (z, w) with |z|, \w\ < 1, and 
if the two branches of the singular variety V meet transversely at (1, 1) as in Figure ^ then for 
some positive constants c and C, 

(L3) a,, = C7 + 0(e-^l('^'^)l) 

for all (r, s) in a certain cone, K, in the positive integer quadrant. This is proved in Theorem K-J.ll 
below, and the constant C computed. Exact statement of the transversality hypothesis requires 
some discussion of the geometry of V. The constant C is computed in terms of some algebraic 
quantities derived from F. The cone K in which this holds is easily described in terms of 
the tangents to the branches of the double pole. The need for some preliminary algebraic and 
geometric analysis to define transversality and to compute C and K motivates our somewhat 
lengthy Section |2j 
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Figure 2. When V has only smooth points 




Figure 3. A pole of order 3 

Example 13.21 below shows the details of this computation for a particular F. The exponential 
bound on the error follows from results in |PemOOj which are cited in Section^ If the multiple 
pole is moved to a point {z*,w*) other than (1, 1), a factor of (z* )"''(«;*)"* is introduced. 

Compare this with the case where V is smooth, intersecting the positive real quadrant as 
shown in Figured In this case one has Ornstein-Zernike (central limit) behavior. Suppose, as 
above, that (1, 1) G V and F has no other poles (z, w) with \z\, \w\ < 1. Then as (r, s) ^ oo with 
r/s fixed, a^s is rapidly decreasing for all but one value of r/s; for that distinguished direction, 
ars ~ C|(r, s)|~^'^ for some constant, C, and the error terms may be developed in a series of 
decreasing powers of |(r, s)\. The two main differences between the double and single pole cases 
are thus the existence of a plateau in the double pole case, and the flatness up to an exponentially 
small correction versus a correction of order |(r, s)|~^/^ for a single pole. D 



Example 1.2 (pole of greater order). Alter the previous example so that V has a pole of some 
order n + 1 at (1, 1), as in FigureEl Then the formula H1.3() becomes instead 

a,, = P(r,s) + 0(e~^l('-'^)l) 

where P is a piecewise polynomial of degree n—1. In Example l3.7l below. P is explicitly computed 
in the case n + 1 = 3. Again, we see the chief differences from the smooth case being a cone of 
non-exponential decay, and an exponentially small correction to the polynomial P{r, s) in the 
interior of the cone. D 

In higher dimensions there are more possible behaviors, but the same sorts of results hold. 
There is a cone on which the exponential rate has a plateau; under some conditions the correc- 
tion terms within this cone are exponentially small. The following table summarizes the results 
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proved in this paper; transversality assumptions have been omitted. 



Theorem number 


Hypotheses 


Asymptotic behavior 


Theorem rm 


two curves in 2-space 


C + exponentially small 


Theorem l'6.'6l 


d + 1 sheets in {d + l)-space 


C + exponentially small 


Theorem 13.61 


more sheets than the dimension 


polynomial + exponentially small 


Theorem 13.91 


fewer sheets than the dimension 


asymptotics start with r "/^"''/^ 


'i'heoreml3.1il 


all sheets tangent 


asymptotics start with r """'/^ 



We conclude the introduction by giving a combinatorial application of Example 11.11 

Example 1.3 (combinatorial application). An independent sequence of random numbers uni- 
form on [0, 1] is used to generate biased coin-flips: if p is the probability of heads then a number 
X < p means heads and x > p means tails. The coins will be biased so that p = 2/3 for the 
first n flips, and p = 1/3 thereafter. A player desires to get r heads and s tails and is allowed to 
choose n. On average, how many choices oi n < r + s will be winning choices? 
The probability that n is a winning choice for the player is precisely 



E 

a+b=n 



n 



(2/3)'*(l/3)' 



r + s 



n 



{l/3Y-^{2/3) 



s-b 



Let Qrs be this expression summed over n. The array {ars}r,s>o is just the convolution of the 
arrays f+') (2/3)'^ (1/3)" and (''+") (1/3)^ (2/3)^ so the generating function F(z, u;) := ^^r 
is the product 



IrsZ^W^ 



F(z,w) 



1 



:i 



1 , 



Applying Theorem 13.11 with G = 1 and detH = —1/9, we see that a^s = 3 plus a correction 
which is exponentially small as r, s — > oo with r/{r + s) staying in any subinterval of (1/3, 2/3). 
A purely combinatorial analysis of the sum may be carried out to yield the leading term, 3, but 
says nothing about the correction terms. The diagonal extraction method of |HK71j yields very 
precise information for r = s but nothing more general in the region 1/3 < r/{r -|- s) < 2/3. D 



2. Preliminary definitions and notation 

For each z G C^^, the truncation (zi, . . . , z^) will be denoted z, and the last coordinate z^+i 
simply by z. We do not specify the size for constant vectors — for example 1 denotes the vector 
(1, . . . , 1)^ of whatever size is appropriate. Thus we write 1 = 1. A minimal singularity z* is 
strictly minimal if V n D(z*) = {z*}. When a minimal point is not strictly minimal, one must 
add (or integrate) contributions from all points of V n T(z*). This step is routine and will be 
carried out in a future article; we streamline the exposition here by assuming strict minimality. 

This article deals entirely with minimal points z* of V near which V decomposes as a union of 
sheets Vj , each of which is a graph of an analytic function z = Uj (z) . The algebraic description 
of this situation is as follows. By the Weierstrass Preparation Theorem /GHTQ there is a neigh- 
borhood of z* in which we may write H{z) = xi'^)W{z) where x is analytic and nonvanishing 
and VF is a Weierstrass polynomial. This means that 



Wiz) 



^n+l 



+ E^^-' 



z)z-' 



j=o 
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where the multiphcity n + 1 is at least 2 if z* is not a smooth point, and the analytic functions 
Xj vanish at z*. 

Now suppose that z* has all coordinates nonzero. Recalling from Lemma 6.1 of |PW02bj 
that any such minimal point of V is locally homogeneous, we see in fact that Xj vanishes to 
homogeneous degree n + 1 — j at z* and that xo has nonvanishing pure z""*" terms for each 
1 < i < d. Then V is locally the union of smooth sheets if the degree n + 1 homogeneous part 
of W (the leading term) factors completely into distinct linear factors, while if V is locally the 
union of smooth sheets then we have a factorization 



(2.1) Hiz)=x{z)l[[z-u,iz)] 

i=o 

for analytic (not necessarily distinct) functions Uj mapping a neighborhood of z* to a neighbor- 
hood of z* . 

Remark. If the expansion of H near z* vanishes to order n + 1 in z, then there are always n + 1 
solutions (counting multiplicity) to H(z,z) = for z near z*. These vary analytically, but may 
be parametrized by n + 1 analytic functions Uj only if there is no monodromy, that is, if one 
stays on the same branch moving around a cycle in the complement of the singular set. Thus a 
multiple point singularity is one whose monodromy group is trivial. We remark also that in two 
variables, local homogeneity of the minimal point z* nearly implies it is a multiple point. The 
only other possibility is a cusp whose tangents are all equal; such a singularity turns out not to 
affect the leading order asymptotics and to affect the lower order asymptotics only in a single 
direction. 

It turns out to be more convenient to deal with the reciprocals Vj = 1/uj. The basic setup 
throughout the rest of this article is as follows. 

Definition 2.1. The point z* of V is a multiple point if there are analytic functions vq, . . . ,Vn,(p 
and a local factorization 



(2.2) F(z) 



n-=o(i-^^.(^))' 

which we call the factored form of F, such that 

• iz*)d+iVj{z*) = 1 for all j; 

• each — ^(z*) / 0; 

• Zd-\-iVj{z) = 1 for some j if and only if z G V. 

Remark. Let Vj denote the local hypersurface parametrized by z = Uj(z). Then the first of 
the conditions says that each Vj passes through z*: there are no extraneous factors in the 
denominator representing surface elements not passing through z* . The last condition says that 
the zeros of the denominator are exactly the poles of F: there are no extraneous factors in the 
denominator vanishing at z* and cancelling a similar divisor in the numerator. 

It may appear that some generality has been lost in imposing the second condition, since we 
are assuming that each sheet of V projects diffeomorphically onto any coordinate hyperplane. 
This latter property is in fact guaranteed by the non- vanishing of the pure z"^ terms of W. 
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To compute (p directly from G and H^ differentiate ()2.1|) n+1 times in the z := z^+i coordinate 
at the point z* to write 

/ ^ X n+l 

We may then write 



-) iJ(z*) = (n + l)!x(z*). 



G G 



ni=o(i - ^^jz) H X n"=o(-^i) n"=o(i - ^v^iz)) 

and solve for (/> at z = z* to obtain 

(n + 1)! G(z*) 



(2.3) ./.(z*) 



*"! J , -, "1" + ! / Q \"+l 



(-(^*)^+l) 



i7(z*) 






For the remaining definitions, fix a strictly minimal element z* € V which is a multiple point, 
and let uq, • • • , Wn be as above. 

Definition 2.2 (Cone of directions corresponding to multiple point). For each sheet Vj and 
multiple point z*, let dirj(z*) be the vector defined by 

dir,(z*):=f^|^,...,^|^,l 

We denote by K(z*) the positive hull of all the dirj(z*) and by Ko(z*) their convex hull, in other 
words the intersection of K(z*) with the hyperplane z^+i = 1- Geometrically, K is precisely 
the collection of outward normal vectors to support hyperplanes of the logarithmic domain of 
convergence of F at the point (log |zj|, . . . , log jz^+il); see |PW02bj for details. 

Let C(z*) be the matrix whose jth row is dirj(z*). We say that z* is nondegenerate if the 
rank of C is d + 1, transverse if the rank is n + 1, and completely nondegenerate if it is 
both transverse and nondegenerate; in this case necessarily n = d and the multiplicity of each 
sheet is 1. 

If iC is a subset of K or Kq consisting of vectors whose directions are bounded away from the 
walls (equivalently the image of K under the natural map ~ : M*^"^^ \ {0} — > MP is a compact 
subset of the interior of K) then we shall say that K is boundedly interior to K (or Kq). □ 

Remark. The importance of K is that analysis near z* will yield asymptotics for r in precisely the 
directions in K(z*). In the smooth case, n = and so K(z*) reduces to a single ray. The point 
z* is transverse if and only if the normals there to the surfaces Vj span a space of dimension n + 1 . 
Alternatively, the n + 1 tangent hyperplanes intersect transversally. When n > d, transversality 
must be violated; nondegeneracy means that there is as little violation as possible. D 

3. Main theorems and illustrative examples 

All of the results in this section are ultimately proved by reducing the problem to the compu- 
tation of asymptotics for a Fourier-Laplace integral (the proofs are presented in later sections). 
Owing to the large number of possibilities arising in this analysis, constructing a complete tax- 
onomy of cases is rather challenging, and the number of potential theorems is enormous. 

We have chosen to present a series of theorems of varying complexity and generality. Taken 
together, they completely describe asymptotics associated with nondegenerate multiple points. 
Our analysis of other types of multiple points requires additional (mild) hypotheses, which will 
almost always be satisfied in applications. The most important case is that of transversal points, 
but we also treat various types of tangencies and degeneracies. However, it is always possible 
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that a practical problem involving a meromorphic generating function may not fit neatly into our 
classification scheme. We hope to convince the reader that in such a situation our basic method 
will yield Fourier-Laplace integrals from which asymptotics can almost certainly be extracted in 
a systematic way. 

In two variables, our results specialize to the following cases. If V is locally the union of 
n + 1 analytic graphs, then either at least two tangents are distinct, in which case Theorem 13.61 
applies, or all tangents coincide. This latter case is more complicated and we require some extra 
hypotheses; see Theorem 13. Ill 

We begin with the simplest case, in which the result admits a relatively self-contained state- 
ment, with as little extra notation as possible. 

Theorem 3.1 (2 curves meeting transversally in 2-space). Let F be a meromorphic function of 
two variables, not singular at the origin, with F{z,w) = G{z,w)/H{z,w) = "^^ ^OrsZ^w^ . 

Suppose that {z*,w*) is a strictly minimal, double point of V. Let 'H.{z*,w*) denote the 
Hessian of H at (z*,w*) and suppose that detii{z* ,w*) ^ 0. 

Then 

(i) for each boundedly interior subset K of'K{z*,w*), there is c> such that 

ars = (^*)-'-K)-^ f , , .,2^^f,2fLr . ,, +0(e-'^'^^'^^')| umformly for {r,s) G K. 

\y/-{z*y{w*)^detll{z*,w*) I 



(ii) if 6 = r/s lies on the boundary of ^.(z* ,w*), then in direction S there is a complete 
asymptotic expansion 

arsr^{z*)-'^{w*)-'Y.^kS-('+'y' 
where bo — '^ 



2-y/-(z*)2(w,*)2detH(2*,«i*)' 



The geometric significance of the nonsingularity of H^z* , w*) is that this is equivalent to the 
curves Vo, Vi intersecting transversally. If, on the other hand, det H(z* ,w*) =0 (equivalently the 
curves are tangent), then higher order information is required in order to compute the relevant 
asymptotics. We treat the latter situation in Theorem 13. Ill 

Note that in the situation of the above theorem, the asymptotic exponential rate of Ors 
is constant on the interior of the cone K(z*,u;*), that is, Urs ~ exp(— (r, s) • v) where v = 
(log z* , log w* ) is constant on K(2:* ,w*); this differs considerably from the asymptotics previously 
derived for smooth points fPW02b _. 

Example 3.2 (a lemniscate). Consider F = l/H, where H{z,w) = 19 — 20z — 2Qw + 5z^ + 
lAzw + bw"^ — 2z'^w — 2zw'^ + z'^w'^. The real points of V are shown in Figure |1J 

The only common zero of H and V H is (1,1), and so this is the only candidate for a double 
point. All other strictly minimal elements of V must be smooth. The second order part of the 
Taylor expansion of F near (1,1) is 4(i(;- 1)^ + 10(z- l)(u;- 1) +4(z - 1)^ = i[w- {3-2 z)][w- 
(3/2 — z/2)]. The Hessian determinant at (1,1) is therefore —36. The two tangent lines have 
slopes —1/2 and —2, so that asymptotics in directions r/s = k G [1/2,2] cannot be obtained 
by smooth points. By Theorem 13.11 for slopes in the interior of this interval, the asymptotic 
Qrs ~ 1/6 holds, whereas on the boundary of the cone (slopes 1/2 or 2), ars ~ 1/12. Directions 
corresponding to slopes outside the interval [1/2, 2] may be treated by the methods of |PWn2bj . 
The exponential order is nonzero in these directions. 
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Figure 4. Figure-eight shape 



This example demonstrates that the local behavior of H at the strictly minimal point (1, 1) 
determines asymptotics. In fact H[l + z,l + w) is homogeneous and irreducible in C[2;, w] and 
hence does not factor globally (in the power series ring C[[z,tt;]]). 

D 

The previous theorem treated the simplest case, when d = n = 1, and our subsequent re- 
sults will be labelled in a similar way. We begin with some theorems involving nondegenerate 
points. The simplest case is when the point is completely nondegenerate (a generalization of the 
hypothesis of Theorem 13. 1|) . 

Theorem 3.3 (d = n; completely nondegenerate). Let F he a function of d + 1 variables with 
F(z) = G{z)/H{z) = ^j.Crz'". Suppose that z* is a strictly minimal, completely nondegenerate 
multiple point ofV and that F is meromorphic in a neighborhood o/T(z*). 

Then there is C such that for each boundedly interior subset K of K(z*), there is c > 
satisfying 

Or ~ (z*)"'" (c + 0(e~'^l''l)^ uniformly for r e K. 
The constant C is given by 

C 



d\i 



,detC(z*)| 

Remark. The picture to keep in mind here is of d-|- 1 complex hypersurfaces in C^"*"^ intersecting 
at a single point. Notice that we again obtain asymptotic constancy on the interior of the cone 
of allowable directions, provided G(z*) 7^ 0. If r belongs to the boundary of K, many different 
asymptotics are possible. D 

In case d < n, a modification of the last result is required, involving some linear algebra which 
we now introduce. Let A = A„ denote the standard n-simplex in M"^-'^, 



ex G 



pn+l 



j=0 



1 and aj > for all j 



> . 



The subspace 1 parallel to A will be denoted by TC. Assuming nondegeneracy of a minimal 
point z* , the rows of C(z*) are n + 1 vectors spanning M "'""^, whence for any 6 G M. ~^^ the space 
of solutions A := A{S) := {a \ aC{z*) = S} will always be (n — d)-dimensional. The set ^ n A 
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is just the set of coefficients of convex combinations of rows of C that yield the given vector S. 
If 5 € Kq then ^ is a subset of the translate of 7i containing A, thus if S is in the interior of 
Kq then ^n A is (n — d) -dimensional. Let A'^ denote the orthogonal complement in Ti of ^O) 
the set A translated to the origin. The dimension (later denoted p) of A'^ is d. The space ^-'" 
does not depend on S, since changing S only translates A. 

The role that A will play is this: there will be an integral over A whose stationary points 
are essentially the points of A, each contributing the same amount; thus the total contribution 
will be the measure of A in the complementary measure which we denote by a. The proper 
normalization in the definition of a is 



a{S) ■.= a{SnAnA) 



fin^diSnAnA) 



/Wn(A) 

where pk denotes fc-dimensional volume in M"+-'^. One may think of this as: a x fi^ = the 
normalized volume measure, /i, on A ( Warning: when ^ is a single point, it is tempting to 
assume that (j{A) = 1, but in fact in this case (j{A) = /i„(A)~^.) 

Proposition 3.4. Assume that the columns of C are linearly independent. Then the projection 
of the column space of C onto A'^ has dimension d. 

Proof. The matrix C is the sum of three columnwise projections, namely one onto Aq, one onto 
A and one onto the span of 1. The sum of the first two projections annihilates the last column 
of C, mapping the column space to a space of dimension d. By definition of A, the space AC is 
a single point, hence AqC = 0, meaning that the first projection is null. Therefore, the second 
projection has rank d. D 

Definition 3.5. When the columns of C are linearly independent, let C denote the matrix 
representing the linear transformation v h^ vC on A'^ with respect to some basis of A'^. Note 
that C is independent of S since A'^ is, and that det C is independent of the choice of basis of 
A'^. For general C, we extend this definition so that C is the projection of the column space of 
C onto the space A , then represented in a (fixed but arbitrary) orthonormal basis of ^ . D 

In stating subsequent results, we shall rely increasingly on derived data such as C. It is 
possible in principle to give formulae for these asymptotics in terms of the original data G and 
H, but such expressions rapidly become too cumbersome to be useful. 

Theorem 3.6 {p = d < n; nondegenerate) . Let F be a function of d + 1 variables, with 
F{z) = G{z)/H{z) = ^j.arz'". Suppose that z* is a strictly minimal, nondegenerate multiple 
point of degree n + 1 ofV and that F is meromorphic in a neighborhood o/T(z*). 

Then K(z*) is a finite union of cones Kj such that for each boundedly interior subset K of 
each Kj, there are c > and a polynomial P of degree at most n — d, such that 

Or = (z*)""" (P(r) + ©(e""^!""!) j uniformly for r e K. 

Here 

^ »(..M^nA) I ^ o((..«)»--) . 

|detC(z*)| 
where A{S) is the solution set to aC{z*) = S with S = vjr^w G Ko(z*), and a is the comple- 
mentary measure. 

Remark. The approximation to a^ is actually piecewise polynomial and is asymptotically valid 
throughout the interior of K. The correction term may, however, fail to be exponentially small 
on some lower-dimensional surfaces in the interior of K where the piecewise polynomial is pieced 
together. 
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We have broken the coordinate symmetry in the formula for the leading term by parametrizing 
r in terms of r^-^-i and S. See section [S] for more comments. D 

Example 3.7 (3 curves in 2-space). The simplest interesting case to which Theorem 13 . 61 applies 
is that where V is the union of 3 curves in 2-space which intersect at a strictly minimal point. 
Suppose that the strictly minimal singularity in question is at (1,1). In this case since d = 1 
and ra = 2, Theorem 13.61 shows that 

ttrs ~ -P(r, s) 

for some piecewise polynomial P of degree at most 1, whenever 5 := r/s lies in the interior of 
the interval formed by the slopes Cj = v'Al). We shall obtain a more explicit expression for P. 

Without loss of generality, we suppose that < cq < ci < C2, at least one of the two 
inequalities being strict. Assume for now that both inequalities are strict: cq < ci < C2- Let 
As = {{o-o, ai, 0:2) G A I oqCq + aici + a2C2 = 6}. When 6 belongs to the convex hull of cq, ci 
and C2, the affine set ^4^ is a line segment whose endpoints are on the boundary of the simplex 
A. Since A is orthogonal to c, the set A is parallel to the projection c of c onto Tl. Letting fi 
denote the mean of the Cj, I!/\/3 the standard deviation, we see that the projection is c — /il 
and its euclidean length is S. 

If Co < S < ci, then one endpoint of the line segment is on the face 02 = 0, with oq = 
-S^,ai = -^^. The other endpoint is when ai = 0, with ao = -^^,02 = -^^- The 
squared euclidean length of this line segment simplifies to 

(Cl -Co)2(c2 - Co)2' 

A similar argument gives the answer when ci < 5 < 02- The squared euclidean length of the line 
segment is then 

(C2 -Ci)2(c2 -Coy 

Both answers agree, and by continuity give the correct answer, when 5 = ci. Since the area of 
A2 is v3/2, the complementary measure a of As is 2/\/3 times its euclidean length. Thus we 
obtain P{r,s) = sf{r/s), where 



Thus 



'(ci-Co)(c2-Co) ' 



if Cq < 5 < ci; 



2? M V, ifci<6<C2. 

, (C2-Cl)(C2-Co) ' ^ — — ^ 



if Co < r/s < ci; 



(3-1) pir,s) = { j'^-:i>tr''^ ' -f ^ / ^ 

l^ (c2^gKc2-co) ' ifci<r/s<C2. 

Note that / = on the boundary of the cone, namely when 5 = cq 01 5 = C2- The above formula 
extends to the case cq = c\ or c\ = C2 in the obvious way. This example illustrates the piecewise 
polynomial nature of the asymptotics in K. The approximation a^s ~ Pi^^ s) is exponentially 
close for b in compact subintervals of (co,ci) U (ci,C2), but requires polynomial correction when 
(^ — > Cl as well as at the endpoints co and C2. □ 

More general formulae for P in this case are described in the forthcoming work |BPn3j . 

We move on to investigate situations which are degenerate according to our terminology. We 
begin with the simplest case, namely when the minimal point is transverse. Note that in this 
case, the set A of linear combinations of the rows of C that yield 8 is always a single point since 
the rows of C are linearly independent. 
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Definition 3.8. If z* is a strictly minimal multiple point and a is an element of A, we define 
the matrix Q = Q{z* , a) to be the Hessian matrix of the function 

e^ -logcxv{z*e'^). 
at the point 6 = 0. We define the matrix M by 

(3.2) M:=M(z%a):=f "^ "^^J^^jV 

y-iC (z ) Q{z ,a)J 

D 

Theorem 3.9 {p = n < d; transverse). Let F be a function of d + 1 variables, with F{z) = 
G{z)/H{z) = ^j.arz'". Suppose that z* is a strictly minimal, transverse multiple point of degree 
n + 1 ofV and that F is meromorphic in a neighborhood o/T(z*). 

Let K be a boundedly interior subset o/Ko(z*) such that detM(z*,Q:(5)) / for S £ K, 
where a{S) is the unique point in A{S). Then there is a complete asymptotic expansion 

Or ~ (z*) '" 2, bk{fd+i)^^ uniformly for S €z K. 

k>0 

If G{z*) 7^ 0, the leading coefficient is given by 

n\ {2TT)^(t>{z*) 



yo 



,/^m detM{z*,a{S)y/^ 
where the square root is the product of principal square roots of the eigenvalues. 

Example 3.10 (2 planes in 3-space). Consider the trivariate sequence {ar^s,t) whose terms are 
zero if any index is negative, and is otherwise given by the boundary condition ao,o,o = 1 and 
the recurrence 

— 2ar^2,s,t — 2ar,s-2,t — (lr,s,t-2- 

Let F{x, y, z) = Ylr s t '^r,s,tx^y^z^ be the associated generating function. Then we obtain 
directly the explicit form 

16 

^"^'^'^^ ~ {4-2x-y-z){4-x-2y-z)' 

The point (1, 1, 1) is a multiple point of V which is readily seen to be strictly minimal. In fact 
this point is on the line segment of multiple points {(1,1,1) + A(— 1,-1,3) : —1/3 < A < 1}. 
We focus here on the point (1, 1, 1) in particular because it is the only point giving asymptotics 
which do not decay exponentially. Near (1, 1, 1), the set V is parametrized by z = 4 — 2x — y and 
z = 4 — X — 2y. The cone K is the positive hull of (2, 1, 1) and (1, 2, 1), with Kq being the line 
segment between these. Other multiple points on the line segment govern other two dimensional 
cones. For example, the point (1/3,1/3,3) corresponding to A = 2/3 governs the cone K2/3 
which is positive hull of (2/3,1/3,3) and (1/3,2/3,3), or in other words, K2/3 is the set of 
(r, s, t) with r + s = t/3 and r,s > t/9. These other cones sweep out, as A varies, a polyhedral 
cone with non-empty interior. There are still directions outside the cone, namely directions 
{r,s,t) with min{r, s} < (r + s)/3; asymptotics in these directions are governed by smooth 
points on one of the two hyperplanes. We return to the examination of the non-exponentially 
decaying asymptotic directions. 
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The multiple point (1,1,1) yields asymptotics for a^st with r + s = 3t and r,s > t. Given 
6 = (2 — a, 1 + a, 1) e Kq, the set A is the single point (1 — a, a). The complement A-^ is always 
all of Ti, which has an orthonormal basis {x}, x := (y^l/2, — -y/l/2). 

The matrix C is equal to ((2, 1) • x, (1, 2) • x) = x. This leads to detM = {A + B + C + D)/2, 
where A, B, C and D are the entries of the restricted Hessian, Q. A routine computation shows 
that A-\-B + C + D = 12, leading to detM = 6 (as is always the case for transverse intersections 
of the pole set, the determinant of M does not vary with direction inside the cone of a fixed 
multiple point). Computing (j){l) from l|2.3j) then gives 



2! 16 

TFT 



^i) = ^^- = i6. 



Plugging this into Theorem 13.91 we find (to first order, as t ^ oo) that 

Or s t ~ ^^^= if r + s = 3f and r/s G (1/2, 2) . 

\/247rt 

The above first order approximation differs from the true value of arst by less than 0.3% when 
(r,s,i) = (90,90,60). D 

As mentioned in the introduction, instead of producing enough variants of these theorems for 
a complete taxonomy, we will stop with just one more result. The case we describe is the most 
degenerate, namely where all the sheets Vj are tangent at z*, so that K(z*) is a single ray. In 
this case det M varies with a and the resulting formula is an integral over a. The methods used 
here may be adapted to prove results when the degeneracy is somewhere between this extreme 
and the nondegenerate cases. 

Theorem 3.11 {p = 0). Let F he a function of d + 1 variables, not singular at the origin, with 
F{z) = G{z)/H{z) = ^j.arz'". Suppose that z* is a strictly minimal, multiple point of degree 
n + 1 ofV and that F is meromorphic in a neighborhood o/T(z*). Further suppose that K(z*) 
is a single ray. 

If detQ{z*,a) / on A, then for r G K(z*) we have 

Or = Z*-'-[6o(r,+l)"-'^/2 ^ o^^^^^^„-d/2-l] _ 

The value of bo is given by 



(27r)rf/2 



/ detQ{z*,a)-^/^dfi{a) 
Ja 



Remark. Theorem 13.111 covers the situation where all the sheets coincide, as opposed to merely 
being tangent. That situation can also be analyzed in a different way by a slight modification 
of the proof of the smooth case, as discussed in _P W02b] . The methods of |PW02bj yield a 
complete asymptotic expansion. D 

Example 3.12. The simplest case illustrating Theorem 13. Ill is that of two curves w = Uj{z) in 
C^, intersecting tangentially at the strictly minimal point z = (1,1). Suppose that Re log Uj(e*^) = 
—dj6'^ + --- with dj > 0. The simplex in question is one-dimensional and we identify it with the 
interval < t < 1. Then for < t < 1 the matrix M(z,f) is just the 1x1 matrix with entry 
df = (1 — t)do + tdi- Hence when do 7^ di we obtain asymptotics in the unique direction 5 of the 
theorem: 
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1,1)^/^ 1 f"' -1/2 , 

y ' dy 



2tt 



\p2M d\ — do y^jj 
2(P{l,l)^Vd^-,/d^ 

\/27r di - do 

0(l,l)Vi 2 



V^ Vdo + Vdi 
This final formula is easily seen to hold also in the case do = di, in which case the formula 
4>{1, l)^/s/^/2^Td^ agrees with the formula in IPWn2bj . D 

Remark. Since analysis in 1 variable is considerably easier than in general, it is possible to derive 
a result for plane curves even when dj vanishes. We omit the details here, but note that such 
a derivation is possible because (by minimality) we must have log wvj{ze) = icj9 — dO"^ + • • • , 
where Cj > 0, m is even and Red > 0. 

If all tangents are equal at {z*^w*)., then 'Vio{z* ,w*) is a singleton. If all ^e\ogwvj{ze^^) 
vanish to the same exact order m and (r, s) G K(z*, ti;*), then 

a,, = {z*)-''-{w*)~'[hos''-^l"' + 0(s"~2/m)] 

with ho given by 



ho 



27r(l - 1/m) di-do 

(t){z,w)V(l/m) 



The above formula for bo also holds in the limit when do = di, when it becomes 6o — , sm 

D 

4. Proofs 

Throughout, we assume that z* is a fixed strictly minimal multiple point of V, of multiplicity 
n + 1, and F{z) = 4'i'^) / YYi=oi^ ~ ■^'^j(^)) near z* as in Definition 12.11 Recall our notational 
conventions s := r^+i, z := z^+i- To ease notation we shall later assume that z* = 1. 

The proofs follow a 4-step process. First we use the residue theorem in one variable to reduce 
the dimension by 1 and to restrict attention to a neighbourhood of z*. Second, the residue sum 
resulting is rewritten in a form more amenable to analysis. Third, the Cauchy-type integral is 
recast in the Fourier-Laplace framework by the usual substitution z = z*e*^. Finally, detailed 
analysis of these Fourier-Laplace integrals (including some generalizations of known results we 
have included separately in |PW02aj to save space here) yields our desired results. 

Restriction to a neighborhood of z*. We first apply the residue theorem in one variable: 
the proof of the following proposition is entirely analogous to that of PW02b, Lemma 4.1]. 

Proposition 4.1 (Local residue formula). Let R{s,z,£) denote the (finite) sum of the residues 
of g : z ^^ z~^~'^F{i, z) at z = Uj(z) inside the ball \z — z*\ < e. 

Let M' be a neighborhood of z* in T(z*). Then for sufficiently small e > 0, 



(4.1) 



ar-{2TTi)-'^ I -z-^-^R{s,z,e)d£ 



is exponentially decreasing in s as s —> oo. 



This is uniform in r as r/s varies over some neighborhood o/Ko(z*). D 
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Rewriting the residue sum. After using the local residue formula H4.1() we must perform 
a d-dimensional integration of the residue sum R. The following lemma, whose proof may be 
found in |DL93l p. 121, Eqs 7.7 & 7.12], will yield a more tractable form for R. Recall from 
Section|21the relevant facts about simplices. In the following, the notation av denotes the scalar 
product X]?=Q OijVj, and h^"-' the nth derivative of h. 

Lemma 4.2. Let h be a function of one complex variable, analytic at 0, and let /x be the 
normalized volume measure on A„. Then 






j=0 



both as formal power series in n + 1 variables vq, . . . ,Vn and in a neighborhood of the origin in 

Corollary 4.3 (Residue sum formula). Let R{s,z,e) be the sum of the residues of the function 
g : z ^^ z^^^^F{z) inside the ball \z — z*\ < e. Define hg^{y) = y*+"0(z; 1/y). 
Then for sufficiently small e, there is 5 > such that for | z — z* | < b, 

(4.2) R{s,z,e)= [ h'f}{a^)dfi{a). 

Proof. First suppose that the functions vq, . . . ,Vn are distinct. Choose e sufficiently small and 
6 > such that g has exactly n + 1 simple poles in \z — z*\ < e whenever |z — z*| < 5. Then 
there are n + 1 residues, the jth one being 



Ur^j (Vr (z) - Vj (z) ) Ur^j (^'r (z) " Vj (z) ) ' 

The result in this case now follows by summing over j and applying Lemma 14.21 

In the case when vq, . . . ,Vn are not distinct, let v*, be functions approaching Vj as t ^ 0, such 
that v*j are distinct for t in a punctured neighborhood of 0, and let g^{z) = z^^^^cl){z)/ YVj=o(^ ~ 
l/v*,(z)), so that g* ^> g as well. The sum of the residues of 5* may be computed by integrating 
5* around the circle |-z — -z*| = £, and since 5* -^ g, this sum approaches the sum of the residues of 
g. Since the expression in (|4.2() is continuous in the variables Vj, this proves the general case. D 

The residue sum formula and Proposition 14.11 combine to show that for some neighborhood 
M' of:?- inT(z*), 
(4.3) 



\{z*Y\ ar-{2iTi)-'^ [ -2-^-1 / /i^"i(Qv)ci/i(a)(^ 
JjV' Ja ^''^ 

where, as usual, we have let s denote r^+i. 



is exponentially decreasing in s, 



Recasting the problem in the Fourier-Laplace framework. In this subsection we assume 
throughout (purely in order to ease notational complexity) that z* = 1. 
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Lemma 4.4 (Reduction to Fourier-Laplace integral). For < k < n, e G M' and ex S M"^^, 
define 

n\{s + n)\ 



^''^^^ • kl{n - k)\{s + k)] 

f(9,a) := - - log(av(e^^)) 
s 



MO, a) ■■= (^) <Pie"',l/y) 



j/=v(e*^) 



I{s;f,i^k) := / e "^li^kd^n+d- 
Then there is a neighborhood J\f of in M , a product of compact intervals, such that 



(4.4) 



(2vr)-'^^^^Pfe(.)/(.;/,V^, 



is exponentially decreasing in s . 



where £ := J\f x A„. 

Remark. Note for later that each -0^ is independent of a, as is /(O, a). 

Proof. Recall that a{A) = nXj^Jn + 1 = Vni^^n)"^ . Applying Leibniz' rule for differentiating h 
yields 

n—k / J \ k 



" / d \^ 



fc=0 

Substituting this into H4.3|) with z = e , shrinking AA' if necessary and letting M be the 
inverse image under 6 h^ e completes the proof and the third step. D 

Analysis of the Fourier-Laplace integraL The proofs of all the theorems in Section |21 
proceed in parallel while we compute the asymptotics for H4.4|) with the aid of the complex 
Fourier-Laplace expansion. To that end, observe that <? is a compact product of simplices 
and intervals of dimension n + d inside M""'" "'"-^, / and ipk are complex analytic functions on a 
neighbourhood of i5 in "H x R'^ (this requires a sufficiently small choice of e in Proposition 14.1(1 , 
and Re f > on £. 

To apply the standard theory, we must locate the critical points of / on £. Recall from 
Definition 12.21 the matrix C = C(z*) whose rows are extreme rays in the cone of outward 
normals to support hyperplanes to the logarithmic domain of convergence of F. 

Proposition 4.5 (Critical points of /). For S := r/s £ Kg, let A{S) denote the solution set 
{a G M"+^ I aC = S}. Let S denote the set of critical points of f on S. Then S consists 
precisely of those points (0, a) with a G A{S) D A. 
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Proof. By strict minimality of z*, for each j the modulus of Vj(z) achieves its maximum only 
when z = z* . Thus any convex combination of Vj (z) with z 7^ z* has modulus less than | Vj (z* ) | . 
In other words, when 6^0, the real part of / is strictly positive, from which it follows that all 
critical points of / are of the form (0, q) for a £ A. 

In fact / is somewhat degenerate: /(O, a.) = for all ct G A, so not only does the real part of 
/ vanish when 6 = 0, but also the Q-gradient of / vanishes there. We compute the 0-derivatives 
at = as follows. For 1 < j < d, 



^^ -^'^' '^-a(^].iz) 



89 j \ s Zd+i \dzj ^ 

Recalling the definition of C, we see that these vanish simultaneously if and only '\i 6 = aC, 
which finishes the proof. D 

The set of critical points of / on if is thus the intersection of £ with an affine subspace of 
]^n+d+i^ We now wish to use known results on asymptotics of Fourier-Laplace integrals. Note 
that in the situation of Theorems 13 . ![ l3^ and l3^ / has a unique critical point on £, whereas in 
Theorem 13.61 S is essentially an afhne subspace of A, and in Theorem 13.1 II tS is essentially just 
A. Though the existing literature on asymptotics of integrals is extensive, we have been unable 
to find in it results applicable to all cases of interest of us here. We have therefore derived the 
extra asymptotics we need in |PWn2aj and cite the relevant ones below in Lemmas 14.71 and 14.81 

Recall that A is the orthogonal complement to (any and every) A{5) in 7i; for each point 
y = (0, a) of S, let By be the product of M with the affine subspace of A containing a and 
parallel to A . 

Then the integrals in Lemma 14.41 mav be decomposed as 

(4.5) / e~''fipkdiin+d= I I e'^ftpkdiip+d] da{y) , 

Js Js \JBy J 

where p := dim^-*-. 

We require a more explicit description of the Hessian matrix of /. 

Proposition 4.6 (Hessian of /). Let y = {0,a) G S. At the point y, the Hessian of the 
restriction of f to By has the block form 



(4.6) M(.'.a)= J_. -''=(^-' . 

\,-iC(z') Q(z-.cc)) 

In this decomposition: 

• the zero block has dimensions p x p, where p = dim^-*- = rankC(z*) 

• the jth column of C{z*) is the projection of the jth column o/C(z*) onto A (expressed 
in some orthonormal basis) 

• Q{z*,a) is the Hessian at y of the restriction of f to the 6-directions, as defined in 
Definition Vl^ 

Proof. Constancy of / in the A directions at 6 = shows that the second partials in those 
directions vanish, giving the upper left block of zeros. Computing (d/d6j)f up to a constant 
gives 

i d ,^ . 

oczj— — V z* 



av(z*) ^ dzj 

and since q;v(z*) is constant when = 0, differentiating in the a directions recovers the blocks 
— iC and —iC . The second partials in the 9 directions are of course unchanged. To see that 



ASYMPTOTICS OF MULTIVARIATE SEQUENCES II 19 



the dimension of A is the rank of C, observe that no nonzero element of A can be orthogonal 
to column space of C, since then it would have been in A; thus the projection of the columns of 
C onto A'^ spans A'^. D 

We may now prove Theorems 13.11 13.31 and 13.91 in reverse order from most general to most 
special. In these cases the function / in the representation (|4.4j) has only one critical point, so 
the decomposition 1)4. 5|) is not really necessary. 

As mentioned above, we have had to derive asymptotic expansions for some cases. We start 
with one of these, applicable to the situations of Theorems 13. H l3^ and 13.91 

Lemma 4.7 (Fourier-Laplace expansion for isolated critical point). Suppose that f has a unique 
critical point £ £, at which the Hessian H is nonsingular. 

(i) // has a neighborhood in £ diffeomorphic to M.'^ then as s ^>- oo there is an explicitly 
computable asymptotic expansion 

oo 
j=0 

In particular 



Oo = -0^(0)- 



7det(H) 

where the square root is defined to be the product of the principal square roots of the 
eigenvalues ofH. 
(ii) // has a neighborhood in £ diffeomorphic to a d-dimensional half-space then there is 
an explicitly computable asymptotic expansion 

oo 

I{s-J,i.,)^Y.^,s-'^"~'"- 

j=0 

Here the value of a^ is half that of the previous case. 

Proof. This is a straightforward corollary of |Hor83l Theorems 7.7.1, 7.7.5, 7.7.17 (iii)]. Details 
are given in |PWn2aj . 

D 

Proof of Theorem 13. 91 We prove the result in the case z* = 1 . The result in the general case 
follows directly via the change of variable z — > z/z* . 

The rank of C is n + 1, so A{d) is a single point and / has a unique critical point, (0,q;). 
Hence A = TC and has dimension n. The critical point is in the interior of £ as long as S is in 
the interior of Kq, which occurs if and only if r is in the interior of K. 

Since .A is a singleton, we may use the formula (|4.4j) and Lemma 14.71 The Fourier-Laplace 
expansion of J^ e"*-* V'fc will have terms s~("-+")/2-i for all j, and when multiplied by the degree 
n — k polynomial pk, will contribute to the s("~'^)/2-'=-i terms for all j. Each of these terms is 
explicitly computable from the partial derivatives of G and H via ipk , but these rapidly become 
messy. The leading term of s("~'^)/2 comes only from A: = 0: 

Or ~ (.^TTj 



y/n + l y^detM{z*;a) 

which proves the theorem in the case z* = 1. D 
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Proof of Theorem 13.31 Here rank(C) = d+ 1 = n+ 1 and so the proof of Theorem 13.91 appHes . 
We can obtam a more exphcit resuh in this case. The block decomposition of M{y) in 1)4. 6() is 
into four blocks of dimensions dxd, implying that M{y) is nonsingular with determinant (det C)^ 
(in performing column operations to bring M{y) to block diagonal form, we pick up a factor of 
(—1)'^ which cancels the factor i'^'^). The last column of C is 1, whose projection onto A'^ is null; 
its projection onto the space spanned by 1 has length \/d + 1, whence | det C| = \/d+ 1 1 det C|. 
Thus we obtain 

Or = C + 0(s"^) 

where 

„ dlcbiz*) 



ldetC(z*)| 

To finish the proof, we quote |FemOOl Corollary 3.2] to see that in fact Or — C is exponentially 
small on boundedly interior subcones of the interior of K. 

Remark. If d is on the boundary of K, it may happen that the point (0, a) has ct lying on a 
face of A of dimension d— 1. In that case, using the second formula from Lemma 14.71 instead of 
the first gives a leading term with exactly half the previous magnitude. Here we do not claim 
exponential decay of Cr — C. In Example 13.21 this occurs when r/s G {1/2, 2}. 

Proof of Theorem 13.11 The second case is covered by the previous remark. 

Details for the first case follow. The analysis of Theorem 13.31 applies. We derive a more 
explicit result, first supposing that {z*,w*) = (1, 1). 

For j € {0, 1}, write Vj{e ) = icjO + 0{9'^), where Cj > 0. The matrix C in this case is simply 
( ^° ^ ). It will be seen below that det C = if and only if det H = 0, so for the purposes of this 
theorem we may and shall suppose that cq / ci. 

We have 

Or-s ~ -. ^ — r uniformly for r/s boundedly interior to [co,ci]. 

|co - Cl| 

This expression for the constant is quite handy, but to recover the statement of the theorem, 
we also express it in terms of the partial derivatives of H. Recall that H = xQ where Q{z, w) = 
(1 — wvo{z)){l — wvii^z)) and x 7^ iisar (1,1). Then at (1,1), Q = Qz = Qw = and so 
Hzz = xQzz = 2xcoCi. Similarly, Hyj^, = x(co + ci) and F^^ = 2x. Thus co,ci are the 
roots of the quadratic {H^u])x'^ — {'^Hwz)x + H^z = 0. Solving via the quadratic formula gives 
|co — ci| = 2\/— detH/H^^. This shows that det H = if and only if cq = ci, as asserted above. 

Formula 1)2. 3() then yields 

_ 0(1,1) _ 2^(1,1) _ G(l,l) 

|co-ci| Hwu,{l,l)\co - ci\ y^-detH(l,l) ' 

The general result follows by a change of scale mapping the multiple point (z* ,w*) to (1, 1). 
Details of the computation are as follows. Make the change of variables z = z/z* ,w = w/w* , 
and write G for G considered as a function of z, w. In an obvious notation, Urs = iz*)'''{w*Yars- 
Then, evaluating at {z,w) = (1,1), we obtain H^ = z*Hz, -ff^g = w*Hw Repeating this yields 
H^ = {z*YH^^, Ha,^ = z*w*H^z and Haa = {w*fHy,yj. 

The special case above now yields 

^ G(l,l) G{z\w*) 

_^^ -^^^ ^/ -{z*Y{w*Y deili{z* ,w*) 

from which the result follows immediately. D 
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To prove Theorems 13 . 61 and 13 . 1 1 1 we require a result on asymptotics of integrals that allows for 
non-isolated critical points. The lemma below says that, as seems intuitively plausible, we may 
compute asymptotic expansions in directions transverse to S and then integrate their leading 
terms along S to get the leading term of the original expansion. 

Lemma 4.8 (Fourier-Laplace expansion for continuum of critical points). Let y £ S and let 
M{y) he the Hessian at y of the restriction of f to By. If ipk{y) 7^ and detM(y) is nonzero 
on S then 

Je \ s J Js Vdet M{y) 

Proof. This is a specialization of jPW02al Thm 2] (the proof of that result is considerably more 
involved than in the isolated critical point case). D 

Lemma 14.81 vields the following expansion: 
(4.7) ar ~ (2vr)-%(.) f ^^ '""'''' / ^M^ daiy). 



s J Js \/det M{y) 

We now proceed to complete the proofs by applying (|4.7j) in each case. We give full details 
only for the case z* = 1; the general case again follows from the change of variable z i— > z/z* 
and we leave the details to the reader. 

Proof of Theorem 13.61 By the nondegeneracy assumption, rank(C) = d + 1, hence p := 
rank(C) = d and dim(By) = 2d. The set S has dimension n — d; we assume this to be strictly 
positive since the case n = d has been dealt with already in Theorem 13.31 At each stationary 
point (0, q) G 5 we find the Hessian determinant of / restricted to By to be the constant value 
(detC)2. Now (jlHl) gives 

ar ~ (27r)-%(s)s-'^(27r)'^,/.(l)ydetAf(y)-i/2rfa(y) 

„_,,/.(l)a(cS) 
~ s 



detCI 



which establishes the formula for the leading term. That the expansion is a polynomial plus an 
exponentially smaller correction again follows from |PemOOl Theorem 3.1]. 

Remark. The expansion via Leibniz' formula is not generally integrable past the leading term, 
and therefore does not give a method of computing the lower terms of the polynomial P. D 

Proof of Theorem I^TTTl Here Kq = {5], and so A{5) = A. Thus dim^-^ = 0, so that for 
y = (0, a) G S, M{y) = Q{z*,cx). Use (|4.7|) once more with dim(By) = d to get 

a. r. i2n)-%{s)s~'/\2nr/' [ -j=4==^Mc^) 

J A y'aetQ{z*,a) 

which simplifies to the conclusion of the theorem. D 

5. Further discussion 

Asymptotics in the gaps. The problem of determining asymptotics when r/s converges 
to the boundary of Kq is dual to the problem raised in jPW02bj of letting r/s converge to 
dJi-o from the outside. Solutions to both of these problems are necessary before we understand 
asymptotics "in the gaps", that is, in any region asymptotic to and containing a direction in the 
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boundary of Kg. For example, in the case Example 13. 121 with (5(1, 1), what are the asymptotics 
for a,^ j,_|_ /p as r ^ oo? 

In addition to the obvious discontinuity in our asymptotic expansions above at the boundary 
of Kq caused by the change from multiple to smooth point regime, in the nondegenerate case 
other discontinuities arise. The leading term may change degree around the boundary of Kq. 
This reflects the behavior of the corresponding Fourier-Laplace integrals, as the affine subspace 
A of positive dimension moves from intersecting the interior of A to only intersecting a lower- 
dimensional face. Subtler effects occur as A passes a corner in the interior of K, when the 
leading term is continuous but lower terms may not be. 

A proper analysis of all these issues in the Fourier-Laplace integral framework requires the 
consideration of integrals whose phase can vary with A. This is the subject of work in progress 
by M. Lladser |L]an3j . 

Effective computation. The greatest obstacle to making all these computations completely 
effective lies in the location of the minimal point z* given r. Assuming the existence of a z*(r) 
with r E K(z*), how may we compute z*(r) and test whether it is a minimal point? Since the 
moduli of the coordinates of z* are involved in the definition of minimality, this is a problem in 
real rather than complex computational geometry and does not appear easy. For example, how 
easily can one prove that the double point in Example 13.21 is minimal? 

Minimal points. Many of our theorems rule out analysis of a minimal point z* if one of its 
coordinates {z*)j is zero. The directions in K(z*) will always have rj = 0, in which case the 
analysis of coefficient asymptotics reduces to a case with one fewer variable. Thus it appears no 
generality is lost. If we are, however, able to solve the previous problem, wherein r/s converges to 
3Ko, then we may choose to let r/s converge to something with a zero component. The problem 
of asymptotics when some rj = o{rk) now makes sense and is not reducible to a previous case. 
Presumably these asymptotics are governed by the minimal point z* still, but it must be sorted 
out which of our results persist when {z*)j = 0. Certainly the geometry near z* has more 
possibilities, since it is easier to be a minimal point (it is easier to maintain \zj\ > \{z*)j\ for z 
near z* when {z*)j = 0). 

In this article and its predecessor we have not treated toral minimal points (minimal points z* 
that are not strictly minimal, and not isolated in T(z*)). Nor have we analysed strictly minimal 
points that are cusps. The extension of our results here to those cases is in fact rather routine, 
and will be carried out in future work. 

Local geometry of critical points. Many of our results fail when the Hessian of / unex- 
pectedly vanishes. Surprisingly, we do not know whether this ever happens in cases of interest. 
Specifically, we do not know of a generating function with nonnegative coefficients, for which 
the Hessian of / on By (notation of Section^ ever vanishes. It goes without saying, therefore, 
that we do not know whether cases ever arise of degeneracies of mixed orders, such as may arise 
in the case of the remark following Example 13.121 

Homological methods. The methods of this paper may be described as reasonably elemen- 
tary though somewhat involved. The necessary complex contour and stationary phase integrals 
are well known, except perhaps for the "resolving lemma" , Lemma 4.2. While writing down these 
results, we found a more sophisticated approach involving some algebraic topology (Stratified 
Morse Theory) and singularity theory. This approach delivers similar results in a more general 
framework. Among other things, it clarifies and generalizes formulae such as the expression (|3.H1 
derived by hand above. Still, the present approach has an advantage other than its chronological 
priority and its independence from theory not well known by practitioners of analytic combina- 
torics. Homological methods, as far as we can tell, are not capable of determining asymptotics 
in "boundary" directions, such as are given in part (ii) of Theorem 3.1. It seems that one needs 
to delve into the fine details of the integral, as we do in the present work, in order to handle 
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these cases. Much of the interest in asymptotic enumeration centers on phase boundaries, so the 
weight of these cases is far from neghgible. We hope some day to integrate the two approaches, 
but at present each does something the other cannot. 

Acknowledgement. We thank Andreas Seeger for giving a reference to an easier proof of 
Lemma 14.21 
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